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Two different controlling methods are proposed to stabi- 
lize unstable continuous-sliding states of a dry-friction oscil- 
lator. Both methods are based on a delayed- feedback mech- 
anism well-known for stabilizing periodic orbits in determin- 
istic chaos. The feedback variable is the elastic deformation. 
The control parameter is either the sliding velocity or the nor- 
mal force. We calculate analytically stability boundaries in 
the space of control parameter and delay time. Furthermore, 
we show that our methods are able to turn stick-slip motion 
into continuous sliding. Controlling friction helps to get a 
better understanding of friction by measuring, e.g., velocity- 
weakening friction forces. 

PACS numbers: 03.20.+i, 46.30. Pa, 07.05.Dz 

If one tries to move two contacting solid bodies lat- 
erally, one often observes stick-slip motion due to dry 
friction (i.e., solid-solid friction with or without lubri- 
cants) M. This motion is characterized by a (more or 
less) periodic switching between sticking (relative sliding 
velocity is zero) and slipping (relative sliding velocity is 
on average much larger than the applied velocity). This 
stick-slip motion is responsible for the everyday expe- 
rience of singing violins and squeaking doors. In most 
technological cases one wants to avoid stick-slip motion 
because it leads to vibrations and wear. The goal is to 
bring the system into the continuous sliding state, where 
the relative sliding velocity is constant and does not os- 
cillate. 
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FIG. 1. The stick-slip oscillator. The feedback variable of 
the proposed controlling methods is the elastic deformation 
a. The control parameter is either the applied velocity v or 
the normal load N. 

The classical method of avoiding stick-slip motion is 
to use a lubricant. But even at high normal load and 
small sliding velocities stick-slip motion occurs M . Here 



we propose an active way to avoid stick-slip motion. It is 
inspired by the methods used in controlling chaos 
The idea is not to change the physics at the friction in- 
terface but to stabilize unstable states. 

In our case the unstable state is the continuous-sliding 
state. It is unstable for velocity-weakening friction laws 
where the friction force decreases with increasing slid- 
ing velocity We use the delayed- feedback method 
proposed by Pyragas for stabilizing periodic orbits in a 
chaotic attractor || . The feedback variable is the elastic 
deformation (see below) . Our parameter of controlling is 
either the sliding velocity or the normal load. 

In order to show that our proposed methods work, we 
have studied analytically as well as numerically a simpli- 
fied model for the lateral motion of two solid bodies in 
contact (see Fig. [I]) . We assume that one of the bodies is 
fixed whereas the other one (with mass M) can slide. The 
elasticity of the sliding bodies (or the whole machinery) 
is modeled by a spring with stiffness k. There is a normal 
load N which presses the bodies against each other and 
there is the applied velocity v which is generally differ- 
ent from the relative sliding velocity x of the bodies. We 
assume that the lateral degree of freedom x is the most 
important one. We therefore neglect all other ones. The 
friction force F at the interface is proportional to N in 
accordance with Amononton's law Q. The equations of 
motion read 

Mx(t) = na(t) - fi K (x(t))N, if x(t) ^ 0, (la) 

or 

i(t)=0, if | /ca(t) | < (isN, (lb) 

and 

a(t) = v - x(t), (lc) 

where /is is static friction coefficient and (j,k(%) is the 
kinetic friction coefficient which in general depends on 
the sliding velocity. The variable a denotes the spring 
elongation (i.e., the difference between the stage position 
and the block position). In general it measures the elas- 
tic deformation of the bodies due to the contact. Our 
feedback variable is a. The control parameter is either 
the applied velocity v or the normal load N. Hence we 
replace either v by 

v = vq + a v [a(t) — alt — r)] (2a) 

or N by 
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N = N + a N [a(t) - a(t - r)], 



(2b) 



where vq and Nq are the unperturbed applied velocity 
and normal load, respectively, r is the delay time, and 
a v and a a? are the amplitudes of control. Note, that load 
control works only if N is always positive. Otherwise it 
would lead to a lift-off of the sliding bodies. We will 
first calculate analytically where in the space of control 
amplitudes and delay time the continuous-sliding state is 
stable. 

The continuous-sliding state is given by x = vq and a = 
Hk(vo)Nq/ k. In order to test its stability we make the 
ansatz x(t) = vq + c± exp(Ai) and a(t) = hk(vq)Nq/ k + 
c a exp(Ai) and linearize the equation of motion in c± and 
c a . Nontrivial solutions are possible only if A fulfills the 
following characteristic "polynomial" : 



MX 2 + [n' K (v )N ~ Ma v (l - e~ Ar )]A + k 
~[^' K MN a v + fJ,K(vo)a N ](l - e 
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The delay term is responsible for e~ Ar which turns the 
polynomial into a transcendental equation for A. Thus 
(|J) has more than two solutions. In fact there are in- 
finitely many. Solutions are either real or coming as con- 
jugated complex pairs. The continuous sliding state is 
stable if the real parts of all solutions of (^) are nega- 
tive. It can be shown that the number of solutions with 
positive real part is finite because the exponential term 
is bounded. In fact one can give upper limits of the real 
part and the imaginary part of the solutions. Actually 
we do not need to known the solutions of (^). We only 
want to know the stability boundary in the space of the 
control amplitudes a v , ajy and the delay time r. The 
stability boundary is a two-dimensional manifold where 
a solution of (j^) crosses the imaginary axis. Such man- 
ifolds can be calculated analytically in parametric form. 
They are solutions of 
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A r = li k (vo)o:n + v' K ( v o)N a v , A^ = Mu;a v , (6) 



and 



B r = k - Muj 2 , B l = n' K {v )N u;. 



(7) 



The parameter ui is the imaginary part of A. Eq. (Q) 
is a linear inhomogencous equation for a v and Qat. By 



giving one a we can express the other one in terms of 
lo. From (||) we get a countable set of solutions for r. 
Near the manifold one can expand A into a Taylor series 
in order to find out whether the real part of A increases 
or decreases when the manifold is crossed. Together with 
the facts that (||) has analytic solutions for a v = = 
and the solutions of (||) are continuous functions of the 
parameters we are able to find the stable regions in the 
parameter space. 
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FIG. 2. Boundaries of stable continuous sliding for (a) pure 
velocity control (i.e., ajv = 0) and (b) pure load control (i.e., 
a v =0). Different boundaries are denoted by the values of 

Figure || shows that unstable continuous sliding states 
can be stabilized with both types of control. Because of 
the multiplicity of the solutions of (Q) we also get stability 
regions for larger values of r. But there are less interest- 
ing because of larger time scales on which the system 
relaxes into the continuous sliding state. Pure velocity 
control (i.e., a at = 0) is less efficient than load control 
because the stability area is smaller and it disappears if 
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FIG. 3. Turning stick-slip motion into continuous sliding. 
The simulations of Eq. (|l|) were done with the purely ve- 
locity-weakening friction law = Mo/(l + x/v a ) with 
fio — 0.5 and v s = 0.5. The simulation starts with an initial 
state near the unstable continuous sliding state. The instabil- 
ity leads to stick-slip motion. At t = 70 the control was switch 
on [velocity control in case (a) and load control in case (b)]. 
The parameters are M = k = (is = No = 1, (a) vq = 0.35, 
a v = —0.33, ajv = 0, r = 4.15, and (b) vo = 0.1, a v = 0, 
a N = -10, r = 0.5. 

In pure load control there is always a stable region. 
Only the minimal value of |cvjv| increases with increasing 
|. For t — > a more simple formula for the stability 
boundary can be given. For small delay times (i.e., r <C 
\J M/k) one can approximate ait) — a(t — r) by the time 
derivative of a. Hence the equation of motion becomes 
a differential equation and (||) becomes a second-order 



pure velocity control (i.e., aN = 0) is not able to stabilize 
unstable continuous sliding states whereas for pure load 
control (i.e., a v = 0) the stabilization works if 



p,' K (v )N 
a N < 7 — ; — • 



(8) 



We have confirmed our analytical results by numerical 
simulations of the equation of motions. Furthermore they 
show that the basin of attraction of a formerly unstable 
continuous sliding state can be quite large. The examples 
in figure || show that stick-slip oscillations disappear after 
the control is switched on. In our simulations we found 
that stick-slip motions survives if the applied velocity v 
is below some critical value v c . Again pure velocity con- 
trol is less robust than load control. For example, stick- 
slip motion can not be destroyed by pure velocity control 
for that value of vo for which in Fig. |3|(b) load control 
turns easily stick-slip motion into continuous sliding. As 
a rule of thumb we found that in the case of pure velocity 
control the sticking time has to be of the same order or 
less than the slipping time. 
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FIG. 4. Boundary of stick-slip motion. The lower 
solid curve is v c for r — > and Coulomb's law (i.e., 
fiK = const < us)- The dotted curves are analytic approx- 
imations (for more details, see the main text). The upper 
solid curve is the lift-off condition. The ratio of both friction 
coefficients is Hk / Us = 0.5. 

In the case of load control it is possible to turn stick- 
slip motion into continuous sliding for arbitrary small 
values of v . To see this we discuss an analytically treat- 
able case. We assume that hk docs not depend on the 
sliding velocity (Coulomb's law). Furthermore we restrict 
ourself to the limit r — > (but cvatt finite). Both assump- 
tions turn the equation of motion into a linear differential 
equation for the slip motion: 



M 



,d x 



dx 

HKCtNT— + k(± - Vq) = 0, (9) 



polynomial since (1 - e ) > At. It is easy to see that ^ thg conditions 
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i(0)=0, ^(0)= M (N + a N Tv ), (10) 

assuming that the slip motion is just starting at t = 0. 
The continuous sliding state is stable for ajv < 0. It will 
be approached in the long-time limit if x(t) > 0, for t > 
(i.e., no resticking). For a^r < —2^/kM/hk, Eq. (||) 
describes an overdamped harmonic oscillator. Thus the 
solution x(t) increases monotonically and never resticks. 
Therefore stick-slip motion disappears even for infinites- 
imal small vq, i.e., v c = 0. In the underdamped case 
one gets resticking if vq < v c . To calculate v c one has 
to solve x(T) = x(T) = 0, with T > 0. This can 
be done numerically. The result for fix — 0.5//g is 
shown in figure ^. In the limits a^r — > and — > 
— 2\fnM I [Ikt we get the approximations v c = v/2^/tF) 
and v c = U7 _1 exp( — 1 — ir^/yl — 7 2 ) , resp., where 
v= (/is — ^k)Nq/\/ kM and 7 = —[IkO-nt kM. Note 
that wo has to be less than Nq/(— a^r) otherwise the 
control mechanism would lead to a lift-off in the sticking 
phase. It is easy to show that the lift-off curve is always 
above v c (see Fig. 

Recently Rozman et al. proposed a different method 
for stabilizing the continuous sliding state Their 
method is similar to the method of Ott et al. for sta- 
bilizing periodic orbits in a chaotic attractor Q . For this 
method one has to reconstruct the Poincare return map 
near the unstable orbit. This is done by observing the 
system dynamics without control. Rozman et al. used 
the normal force as the parameter of controlling. 

Compared with our method the advantage of the 
method of Rozman et al. is that one has not to rely 
on macroscopic equations of motion like ([!]) . Such equa- 
tions of motion are reliable for large sliding velocities but 
it is well-known that they may be not correct for small 
velocities, especially in the case of transitions from stick- 
ing to sliding and vice versa |.| [l(J. In fact Rozman et 
al. have tested their method for a simple model where 
in addition to the macroscopic degree of freedom (i.e., 
the position of the sliding block) an internal degree of 
freedom appears which describes the state of a lubricant. 

There are two disadvantages of the method of Roz- 
man et al.. First, controlling methods a la Ott et al. 
work only in the vicinity of periodic orbits. If these or- 
bits are embedded in a chaotic attractor the system will 
eventually come close to them. Therefore, turning stick- 
slip motion into continuous sliding is possible only, if the 
stick-slip motion is erratic enough to be close to the con- 
tinuous sliding state. Otherwise, the method works only 
if one starts at a large stage velocity where the continuous 
sliding state is already stable and then slowly decreases 
the velocity below the value where the continuous sliding 
state becomes unstable Q. A delayed- feedback method 
is not restricted to the vicinity of the unstable orbit. Of 
course starting far away from the orbit may lead to strong 
controlling forces at the beginning (see Fig. |). But they 



decay exponentially by approaching the orbit. The sec- 
ond disadvantage of the method of Rozman et al. is the 
necessity of reconstructing the dynamics. This may be 
more or less difficult depending on the details of the dy- 
namics of the internal degrees of freedom at the friction 
interface. For technological applications this might be 
important especially because the reconstruction has to 
be recalibrated from time to time. 

In this Letter we have introduced two robust methods 
of stabilizing continuous sliding. They are also able to 
destroy regular stick-slip motion. Both methods rely on 
a delayed feedback where the feedback variable is the 
elastic deformation of the sliding bodies or the machinery. 
The controlling parameter is either the applied velocity 
or the normal load. The velocity control is less robust 
than the load control. 

There are two fields of application of controlling fric- 
tion. Obviously there will be technological applications 
for reducing vibration and wear. But controlling fric- 
tion experiments can also be used to increase our un- 
derstanding of the physics of dry friction. For example, 
using these methods one can measure the effective fric- 
tion force as a function of the sliding velocity even in the 
velocity-weakening regime. 
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